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I review the status of theoretical aspects of B-decays. The principal difficulty in interpreting the wealth of 
experimental data is the control of non-perturbative QCD effects, and the talk is focused on attempts to 
control these effects. Lattice results for the decay constants, B-B mixing and semileptonic form-factors are 
summarized. The discrepancy of the theoretical predictions and experimental measurements for the ratio of 
lifetimes t(A;,)/t(_Bo) is discussed, as well as the status of the semileptonic branching ratio of the B-meson. 
The difficulties in making quantitative predictions for exclusive nonleptonic decays are stressed, and some 
recent approaches to this problem are outlined. 



1 Introduction 



Weak decays of heavy quarks are a particularly fertile field of research for detailed tests of the 
Standard Model of Particle Physics, for measurements of its parameters (the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements in particular) and for potential signatures of new physics. During 
the talk of P. DrellEl we have seen many exciting new experimental results on B-decays, and the flow of 
new data will continue for many years from existing and new facilities. The main theoretical difficulty 
in interpreting the experimental data is the control of the non-perturbative strong interaction effects, 
and this problem is the main focus of my lecture. 

For some of the physical quantities discussed below I will summarize status of Lattice results; 
these results have been taken from the recent review written together with J. FlynnEI, where more 
details can be found. The lattice formulation of QCD, together with large scale numerical simula- 
tions, enables one to calculate non-perturbative QCD effects, from first principles, with no model 
parameters or assumptions. The precision of the calculations is limited, however, by the available 
computing resources, and much effort is being devoted to reducing the systematic errors. These 
uncertainties, and the theoretical and computational efforts to control and reduce them, have been 
discussed at this conference in some detail by M. LuscherB. I will therefore limit my comments to 
the specific computations I am discussing. 

The plan for the talk is as follows. In the next section I will review the status of lattice 
calculations of the leptonic decay constants of heavy mesons and of the amplitudes for the important 
process of B-B mixing. Sec. [| contains a review of exclusive semileptonic decays of B-mesons which 
are being used to extract the V u b and V c b CKM matrix elements. Inclusive semileptonic decays 
are discussed very briefly in sec. |], where the recent suggestions to determine V u b by measuring 
the hadronic invariant mass spectrum are outlined. I then digress from the mainstream of the 
presentation, to discuss the problem of power corrections to hard scattering and decay processes in 
general, and in the heavy quark effective theory (hqet) in particular (sec. ^|). The next two sections 
contain discussions of inclusive (sec. ||) and exclusive (sec. ^) nonleptonic decays. In particular in 
sec. |^ I will discuss two topics which have received much attention lately, viz. the lifetimes of beauty 
hadrons and the semileptonic branching ratio of B-mesons. Finally sec. ^contains some conclusions. 
CP-violation in B-decays, which is perhaps the principal motivation for the coming generation of 
^-factories, is the subject of the talk of Y. NirQ, and will not be discussed below. 
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2 Leptonic Decays and B a —B° Mixing 

I start with a discussion of the calculations of the leptonic decay constants and of the B-parameter 
of B-B mixing. 

2.1 Leptonic Decay Constants 

The decay constant of a meson is a single number which contains all the information about the 
non-perturbative QCD effects in leptonic decays of the meson (see e.g. Fig. [l]). Parity symmetry 
implies that only the axial component of the V — A weak current contributes to the decay, and 
Lorentz invariance that the matrix element of the axial current is proportional to the momentum of 
the meson (with the constant of proportionality defined to be the decay constant): 

(0|4,(0)|fl(p)>=t/ B jv, (1) 

and similarly for the £>-meson. Knowledge of fs would allow us not only to predict the rates for 
leptonic decays, but would also be very useful for describing B-B mixing as explained below, as well 
as for our understanding of other processes in £?-physics, particularly those for which "factorization" 
turns out to be a useful approximation. 

b 



B~ 

u 

Figure 1: Diagram representing the leptonic decay of the B- meson. 

Many lattice groups have evaluated, or are evaluating, fjj and fg. Our view of the current 
status of the calculations can be summarized by the following values for the decay constants and 

tllOll* I*El1j10S Q* 

fn = 200±30MeV, f Ds = 220 ± 30 MeV , 

f B = 170±35MeV, f Bs — 195 ± 35 MeV , (2) 

fojfo = 1.10 ±0.06, f B Jf B = 1.14 ±0.08 . 

All the results are presented using a normalization for which ~ 131 MeV. 

In the absence of experimental results on the decay constants of the S-meson, the results above 
represent our best estimates of this quantity. During the past three years or so, experimental mea- 
surements of fjj B have been made by several groups, and at this years Heavy Flavours conference 
the raporteur summarized the results by: 

fn, = (251 ± 30) MeV . (3) 

The agreement of the lattice predictions with the experimental result is very satisfying, and gives us 
confidence in the results for the other decay constants. 

The results for the decay constants have generally remained stable for many years, but the 
quoted errors have not decreased significantly in this time. This is because the errors are dominated 
by systematic effects, and it is difficult to decrease the uncertainties due to the quenched approxima- 
tion (i.e. the neglect of vacuum polarization effects) without performing reliable calculations with 
dynamical fermions, which will still take several more years. For example, the value of the lattice 
spacing typically varies by 10% or so depending on which physical quantity is used to calibrate the 
lattice simulation. This variation is largely due to the use of the quenched approximation. We there- 
fore consider that ±10% is an irreducible minimum error in decay constants computed in quenched 
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Figure 2: Continuum extrapolation of JLQCD U results (preliminary) for fo and fs- The graph shows agreement 
between results from two different lattice formulations when extrapolated to zero lattice spacing. Open symbols 
denote an unimproved Wilson action (c sw = 0), filled symbols denote an improved Shcikholcslami-Wohlert or Clover 
(c S w = 1) action. 

simulations (and a correspondingly larger error when physical quantities of higher dimension are 
calculated). The remainder of each error in_eqs. (|J) and in the lattice results quoted below is based 
on our estimates of the other uncertainties!!], particularly those due to discretization errors and the 
normalization of the lattice composite operators. Much successful work continues to be done to re- 
duce these uncertainties □, and the results have remained stable during these improvements, adding 
significantly to our confidence, but not reducing the quoted errors by very much. 

As an example of the efforts of some recent work aimed at reducing the systematic uncertainties 
in the evaluation of the decay constants consider Figure [2] which contains preliminary results from 
the JLQCD collaborations. The figure shows results obtained at several different values of the lattice 
spacing a and using two lattice formulations of the heavy quark action, the standard Wilson for- 
mulation (but with a modified normalization of the fields, using the so-called KLM-normalizationQ) 
and with the SWB or "clover" improved action. The lattice was calibrated using the string tension 
to determine the spacing. The feasibility of obtaining (statistically) accurate results for a series of 
lattice spacings and actions is a new development, largely made possible by the increase in com- 
puting power. The figure also shows extrapolations to the continuum (a — 0) limit. In this case 
the continuum values of fs and fo obtained using the two formulations agree remarkably well; the 
agreement is still acceptable (although not so remarkable) when quantities other than the string 
tension are used to determine the lattice spacing. We also note that in this case the dependence 
on the lattice spacing is milder for the improved action as expected (although further studies are 
needed to confirm whether this is a general feature). 

The MILC collaboration has began to study the effects of quenching, and their very early findings 
are that the decay constants may be a little larger (perhaps by about 10%) when the effects of quark 
loops are included □. Further work will be done in the next few years to study this further. 

2.2 Bb -Parameter 

In order to obtain information about the unitarity triangle from experimental studies of B-B mixing, 
one needs to control the non-perturbative QCD effects which are contained in the matrix element 
of the AB=2 operator: 



M(/i) = (B°\ 67m(1-75)9 &7 M (l-7s)<? \B°), 



(4) 
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where the light quark q=d or s The argument /i implies that the operator has been renormalized at 
the scale p.. It is conventional to introduce the i?s-parametcr through the definition 

M(n) = ^f%m%B B {n). (5) 

The dimensionless quantity Bb is better-determined than f B in lattice calculations, so that the 
theoretical uncertainties in the value of the matrix element M , needed for phenomenology, are 
dominated by our ignorance of 

Bb (p) is a scale dependent quantity for which lattice results are most often quoted after trans- 
lation to the MS scheme. The next-to-leading order (NLO) renormalization group invariant B— 
parameter (5^°) is defined by 



6$° = a s ( M )- 2 /* (l + ^-Jn f ] BbM, (6) 

where (3q = 11 — 2rj//3_and J nj is obtained from the one- and two-loop anomalous dimensions of 
the AB=2 operator byEJ, 

j -^( 70 !^ 71 )' (7) 

with {3\ — 102 — 38rt//3, 70 = 4 and 71 = —7 + 4n//9. In the discussion below we choose /i = to/,. 

A number of groups have evaluated B in quenched lattice simulations E , from which we deduce 
the preferred value □ 

Bt = 1.4(1) • (8) 
The relevant quantity for B-B mixing is f B Bs- Taking the result in eq. (^) above for B B l ° with 
f B = 170(35) MeV from eq. (0) gives 



fs = 201 (42) MeV (9) 

as our lattice estimate. For the phenomenologically important quantity £, which relates the ampli- 
tudes for Bd — Bd and B s — B s mixing, we findH: 



$ _ !»J&L = 1 .14(8). (10) 
/fid v B Bd 

3 Exclusive Semileptonic Decays of B-Mesons 

Semileptonic decays of B-mesons are currently being used to determine the V c b and V u b CKM-matrix 
elements. Exclusive decays are represented by the diagram in Figure ||. It is convenient to use space- 
time symmetries to express the matrix elements in terms of invariant form factors (using the hclicity 
basis for these as defined below). When the final state is a pseudoscalar meson P, parity implies 
that only the vector component of the V — A weak current contributes to the decay, and there are 
two independent form factors, / + and /°, defined by 



(P(k)\V»\B(p)) = f + (q 2 ) 



r 



/V) B 2 P Q», (11) 



where q is the momentum transfer, q = p—k. When the final-state hadron is a vector meson V, 
there are four independent form factors: 

(V(k,e)\V»\B(p)) = 2V{ f ] e^e^h (12) 



(V(k, e)\A»\B{p)) = i(m B +m v )A 1 (q 2 )e* " - i A ^ Q ) e*-p {p+ky + i^-l2m V £*-pq»,(13) 

rriB+mv q A 



4 



where e is the polarization vector of the final-state meson, and q = p—k. Below we shall also discuss 
the form factor Aq , which is given in terms of those defined above by Aq = A + A3 , with 




<1 

Figure 3: Diagram representing the semileptonic decay of the _B-meson. q represents the light valence antiquark, and 
the black circle represents the insertion of the V— A current with the appropriate flavour quantum numbers. 



3.1 B -> D* and B -> D Decays 

Exclusive B — > D* and, more recently, B — > D semileptonic decays are used to determine the V c b 
element of the CKM matrix. Theoretically they are relatively simple to consider, since heavy quark 
symmetry implies that the six form factors are related, and that there is only one independent form 
factor specifically: 



f + (q 2 ) = V(q 2 ) = A (q 2 ) = A 2 (q 2 ) 



1 



Q 



(M B + M D ) 2 



Mr 



Mr 



2sJM B M D 



(15) 



where to = vb • vd- Here the label D represents the D- or _D*-meson as appropriate. In this leading 
approximation the pseudoscalar and vector mesons are degenerate. The unique form factor £(w), 
which contains all the non-perturbative QCD effects, is called the Isgur-Wise (IW) function. Vector 
current conservation implies that the IW-function is normalized at zero recoil, i.e. that £(1) = 1. 
This property is particularly important in the extraction of the V c b matrix element. 

The relations in eq. ( |l5| ) are valid up to perturbative and power corrections. The precision 
with which V c b can be extracted is limited by the theoretical uncertainties in the evaluation of these 
corrections. Nevertheless we are in the fortunate situation that it is uncertainties in corrections 
(which are therefore small) which control the precision. 

The decay distribution for B — ► D* decays can be written as: 



dr 



G, 



duj 48tt 3 



{M B -M D ,fMl, x/uj 2 - 1 [uo+lf 



4cu M% - 2wM B M D . + M 2 D , 2 2 

1 + ^TT {Mb — Mb*) 2 J ^ T M ' 

(16) 

where T(iS) is the IW-function combined with perturbative and power corrections. It is convenient 
theoretically to consider this distribution near u = 1. In this case ^(1 ) = 1, and there are no 
Oil/mo) corrections (where Q — b or c) by virtue of Luke's theorem E3, so that the expansion of 
^"(1) begins like: 



T(l)=VA 1 + 0- 



C2 



A 2 

QCD 



(17) 



where tja represents the perturbative corrections. The one-loop contribution to t\a has been known 
for some time now, whereas the two- loop contribution was evaluated last year, with the result E3: 



riA = 0.960 ±0.007 



(18) 
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where we have taken the value of the two loop contribution as an estimate of the error. 

The power corrections, .are much more difficult to estimate reliably. Neubert has recently com- 
bined the results of refs. oli3ii3 to estimate that the Oil /nig) terms in the parenthesesin eq. (|l7j ) 
are about -0.055 ± 0.025 and that = 0.91(3) 0. Bigi, Shifman and UraltsevB, consider 

the uncertainties to be bigger and obtain 0.9 1(6). —Combining the latter, more cautious, theo- 
retical value of !F(1), with the experimental resultta .F(l)|T4f,| = (34.3 ± 1.6)10~ 3 , readily gives 
\V c b\ = (37.7 ± l.&exp ± 2.5 t ft)10" 3 . In considering this result, the fundamental question that has to 
be asked is whether the power corrections are sufficiently under control. This will be discussed in 
more detail in section |^. 

Having discussed the theoretical status of the normalization let us now consider the shape 

of the function T{uj), near uj — 1. A theoretical understanding of the shape would be useful to guide 
the extrapolation of the experimental data, and also as a test of our understanding of the QCD 
effects. We expand T as a power series in u> — 1: 

T{u) = JF(1) [l - p 2 (w - 1) + c (w - l) 2 + ■ ■ •] , (19) 

where 

P 2 = P 2 + (0-16 ± 0.02) + power corrections , (20) 

and p 2 is the slope of the IW-function. 

Recently, using analyticity and unitarity properties of the amplitudes, as well as the heavy 
quaxk symmetry Caprini and Neubert have derived an intriguing result for the curvature parameter 
c EHl, c ~ 0.66 p 2 — 0.11. This result would have effectively removed one of the parameters from the 
extrapolation of the experimental data to u> = 1. The derivation of this relation has been criticized 
in ref. c3, primarily for discarding sub-threshold (B c ) contributions, and the expectation is that 
the corrected relations give somewhat weaker information. Caprini et al_. are preparing a revised 



slope-curvature relation Efl; preliminary indications were presented in 

Finally in this section I consider B — > D sepileptonic decays, which are beginning to be measured 
experimentally with good precision (see refs.l^il'cJ and references therein). The decay distribution 
can be written as: 

- ^{Mb + m d YmUJ - i) f |v; h | 2 s 2 M , (2i) 

where again Q{u>) is the IW-function combined with perturbative and power corrections. Theoreti- 



cally the first complication is t 
out by Shifman and Voloshin 



t the 1/rriQ corrections to (7(1) do not vanish. However, as pointed 
, these corrections would vanish in the limit in which the b and 
c-quarks are degenerate, and hence are suppressed. Ligeti^Nir, and Neubert estimate the 1/mq 
corrections to be between approximately —1.5% to +7.5%I 2 J. The 1/toq corrections for this < 
have not yet been studied systematically. A recent theoretical estimate for C/(l) gives 0.98 ± 0.07 

3.2 B — > p and B — > it Decays 

In this subsection we consider the heavy-to-light semileptonic decaySj-Bj^ p and B — > ix which are 
now being used experimentally to determine the V u b matrix element l^a. Heavy quark symmetry 
is less predictive for heavy-to-light decays than for heavy-to-heavy ones. In particular, there is no 
normalization condition at zero recoil corresponding to the relation £(1) = 1, which is so useful in the 
extraction of V c b- The lack of such a condition puts a premium on the results from nonperturbative 
calculational techniques, such as lattice QCD or light-cone sum rules. Heavy quark symmetry does, 
however, give useful scaling laws for the behaviour of the form factors with the mass of the heavy 
quark (toq) at fixed u>: 

f+,A Q ,A 2 ,V ~ ^h-Q- ; A 1 ,f°^^ = ; A 3 ~ m Q . (22) 
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These scaling relations are particularly useful in lattice simulations, where the masses of the quarks 
are varied. Moreover, the heavy quark spin symmetry relates the B —> V matrix elements Ej'EH 
(where V is a light vector particle) of the weak current and magnetic moment operators, thereby 
relating the amplitudes for the two physical processes B° — > p + l~vi and B — > K*j, up to SU(3) 
flavour symmetry breaking effects. These relations also provide important checks on theoretical, and 
in particular on lattice, calculations. 

Recent work includes detailed lattice studies by several groups, but in particular by the UKQCD 
collaboration, who try to exploit all possible symmetries in order to optimize the available informa- 
tion about the form factors and calculations using light-cone sum rules. 

3.3 Lattice Calculations of B p,ir Semileptonic Decays and of the Rare Decay B — > K*j 

The techniques required to extract the form- factors for B — > ir + l~vi and B — ► p + l~i>i semileptonic 
decays are very similar to those used to compute the short distance contribution to the rare radiative 
decay B — > K*j, so I consider them together. For completeness, I define here form factors for the 
matrix element of the magnetic moment operator responsible for this decay: 

3 

(K*(k,e)\sa^b R \B(p))=Y,Ci 1 T t (q 2 ), (23) 

where q = p—k, e is the polarization vector of the K* and 

C\ = 2e fw > p e*»p x k<>, (24) 
Cl = el(m 2 B -m 2 K ,)-e-q(p + k)^ (25) 

Cl = e*.qU- g ' (p + k)X (26) 

V rn B m K* / 

T3 does not contribute to the physical B — > K*-y amplitude for which q 2 = 0, and Ti(0) and T2(0) 
are related by, 

T 1 {q 2 =Q)=iT 2 (q 2 =Q). (27) 

Hence, for the process B —> K*-f, we need to determine T\ and/or T2 at the on-shell point q 2 =0. 

From lattice simulations we can only obtain the form factors for part of the physical phase space 
for all the decays. In order to control discretization errors we require that the three- momenta of the 
B, 7r and p mesons be small in lattice units. This implies that we determine the form factors at large 
values of momentum transfer q 2 = (job — P-k, p ,k*) 2 - Experiments can already reconstruct exclusive 
semileptonic b — ► u decays (see, for example, the review inE3) and as techniques improve and new 
facilities begin operation, we can expect to be able to compare the lattice form factor calculations 
directly with experimental data at large q 2 . A proposal in this direction was made by UKQCD EJ 
for B° — > p + l~vi decays. To get some idea of the precision that might be reached, they parametrize 
the differential decay rate distribution near q^^ by: 

?0 _^ n +l-r,.\ r$ IT/ . |2 



dT(B° -> p+l-Dj) = 12 G 2 F \V ub 
dq 2 192tt 3 M| 



10-12 ToStFT I" Ai (? 2 ) « 2 (1 + % 2 -9LJ) , (28) 



where a and b are parameters, and the phase-space factor A is given by X(q ) = (m B + m p — q ) — 
Am 2 B m 2 . The constant a plays the role of the IW function evaluated at u = 1 for heavy-to-heavy 
transitions, but in this case there is no Sianmetry to determine its value at leading order in the heavy 
quark effective theory. UKQCD obtained 

a = 4.6+8J±0.6GeV and b = (-8+ 1) x lO" 2 GeV 2 . (29) 
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q 2 (GeV 2 ) 

Figure 4: Differential decay rate as a function of q 2 for the semileptonic decay B° — > p + Z~P;, taken from0. Points 
are measured lattice data, solid curve is fit from eq. with parameters given in eq. (|2q|) . The dashed curves show 
the variation from the statistical errors in the fit parameters. The vertical dotted line marks the charm endpoint. 



The fits are less sensitive to b, so it is less well-determined. The result for a incorporates a systematic 
error dominated by the uncertainty ascribed to discretization errors and would lead to an extraction 
of \V u b\ with less than 10% statistical error and about 12% systematic error from the theoretical 
input. The prediction for the dT/dq 2 distribution based on these numbers is presented in Fig. ||. 
With sufficient experimental data an accurate lattice result at a single value of q 2 would be sufficient 
to fix \V ub \. 

In principle, a similar analysis could be applied to the decay B° — > tt + 1~vi. However, UKQCD 
find that the difficulty of performing the chiral extrapolation to a realistically light pion from the 
unphysical pions used in the simulations makes the results less certain. The B — > tt decay also has a 
smaller fraction of events at high q 2 , so it will be more difficult experimentally to extract sufficient 
data in this region for a detailed comparison. 

We would also like to know the full q 2 dependence of the form factors, which involves a large 
extrapolation in q 2 from the high values where lattice calculations produce results, down to q 2 = 0. 
In particular the radiative decay B —> K*j occurs at q 2 — 0, so that existing lattice simulations 
cannot make a direct calculation of the necessary form factors. In order to determine the form 
factors at lower values of q 2 from their measurements the lattice collaborations, and the UKQCD 
collaboration in particular, exploit a number of important constraints. I now briefly outline these 
in turn: 

• An interesting contribution to the problem of extrapolation to low q 2 has been suggested by 
LellouchEj for B° — > ir + l~i)i decays. Using dispersion relations constrained by UKQCD lattice 
results at large values of q 2 and kinematical constraints at q 2 =0, one can tighten the bounds 
on form factors at all values of q 2 . This technique relies on perturbative QCD in evaluating 
one side of the dispersion relations, together with general properties of field theory, such as 
unitarity, analyticity and crossing. It provides_model-independent results which are illustrated 
in Fig. |. The results (at 50% CL - see Ref.E3 for details) are 

/+(0) = 0.10-0.57 and T{B° -> tt+Z"^) = 4.4-13 \V 2 b \ ps" 1 . (30) 

Unfortunately these bounds are not very restrictive when constrained by existing lattice data. 
In principle, this method can be applied to B — > p decays also, but is, more complicated there, 
and the calculations have yet to be performed. Recently, Becirevicc3 has applied the method 
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q 2 (GeV 2 ) 



Figure 5: Bounds on /"*" and / for B — > ir + l ui from dispersive constraintsEj. The data points are from UKQCDEJ, 
with added systematic errors. The pairs of firre-icurves are, outermost to innermost, 95%, 70% and 30% bounds. The 
upper and lower shaded curves are light-cone a and three-point EJ sum rule results respectively. 



for the process B — > K*j, using lattice results from the APE collaboration as constraints. 
However, he has not applied the kinematic constraint at q 2 — and the resulting bounds are 
not informative: they become so, however, once he uses a light-cone sum rule evaluation of the 
form-factors for the process B — > K*j as an additional constraint. These dispersive methods 
can be used with other approaches in addition to lattice results and sum rules, such as quark 
models, or even in direct comparisons with experimental data, to check for compatibility with 
QCD and to extend the range of results. 

• UKQCD make use of the kinematic constraints on the form factors at q 2 — 0: 

/+(o) = /°(o), r 1 (o) = iT 2 (o), 4,(0) = 4,(0). (3i) 



In spite of all the constraints, model input is required to guide q 2 extrapolations. We can 
ensure that any assumed q 2 -dependence of the form factors is consistent with the requirements 
imposed by heavy quark symmetry, as shown in (|22] ) , together with the kinematical relations of 
eq. (|3l]). UKQCD also verify that the expected relations between form factors of the different 
processes (at fixed w) 

4(<Z 2 M) - 2iT 2 (g a (u;)), V(q 2 (oj)) = 2T 1 (q 2 (oj)), (32) 

are indeed well satisfied in the infinite mass limit. Even with all these constraints, however, 
current lattice data do not by themselves distinguish a preferml_q 2 -dependence. Fortunately, 
more guidance is available from light-cone sum rule analyses E3'L3 which lead to scaling laws 
for the form factors at fixed (low) q 2 rather than at fixed Li as in eq. (|2^). In particular all 
form factors scale like M~ 3 / 2 at q 2 = 0: 



(33) 



where / labels the form factor, M is the mass of the heavy-light meson and is a calculable 
leading logarithmic correction. 
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A combined fit to the lattice data for the form factors (obtained at large values of q 2 ) satisfying all 
these constraints is shown in Fig. The quark masses have been chosen to correspond to the K* 
vector-meson. The figure demonstrates the large extrapolation needed to reach q 2 = 0. 

_Dur preferred results for B° — > ir + l~Di and B° — > p + l~T>i come from the UKQCD constrained 
fitsEa. Their best estimates of the total rates are: 

T(B° ^7r+e-D e ) =8.5±f:||V^ 6 | 2 pa- 1 and T{B° -> P +e-D e ) = 16.5+ 1;§ \V ub \ V" 1 • (34) 



There are also preliminary results for heavy-to-light, form factors from FNAL, JLQCD and a 
Hiroshima-KEK group (see the reviews inEJ'El) and the different lattice calculations are in agreement 



for hra? 
n@0) 

for the form factors at large q 2 where they are measured. 

For the form factor Ti(0) of B — > K*-f decays, the combined fits give 



R «* = z w^w = ^) |T(0)|2 



Ti(0) = 0.16(f). (35) 

Using—this value to evaluate the ratio (given at leading order in QCD and up to 0(1 /m 2 ) correc- 
tions!^) 

r(S^if* 7 ) Jm B 

rt K * = 

results in 

R K , = 16(f)%, (37) 

which is consistent with the experimental result 18(7)% from CLEoEl. Discrepancies between Rk' 
calculated using T(0) and the experimental ratio T(B — > K*j)/T(b — > 57) could reveal the existence 
of long-distance effects in the exclusive-decay. It has been proposed that these effects may be 
significant for the process B — > if*7E3 ~e3. but within the precision of the experimental and lattice 
results, there is no evidence for them. 
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Figure 7: Semileptonic B form-factors as functions of t = q 2 obtained from light-cone sum rules ~ . The dashed lines 
give error estimates and the points are from the lattice simulation of the UKQCD collaboration 



3.4 Light Cone Sum Rule Studies of the Form Factors for B° — ► vi and B — > p + l vi decays 

The second-technique which has been used to study heavy —> light decays is QCD sum rules. Ball 
and BrauncH have recently clarified the origin of the discrepancy iaxesults obtained using different 
types of sum rules for B — > p semileptonic decays at large recoil e3'c3. These authors explain why 
the standard sum rules, in which one performs an Operator Product Expansion (OPE) in terms of 
operators of increasing dimension, fail in this region of phase space. They stress the necessity of 
using light cone sum rules, in which the contributions are classified by the components of different 
"twist" of the distribution amplitude xtfjthc p-meson (strictly speaking of its moments) . This work 
exploits and extends earlier studiest3E§CJ. A careful investigation of the dominant contributions to 
the form factors at large recoil yields the scaling law in eq. ( |33| ) for the behaviour of the form factors 
with the mass of the heavy quark at small fixed q 2 (rather than fixed u>). 

The form factors derived in ref. ED are shown in Fig. (7] together with the results from the 
UKQCD collaboration at large values of q 2 . In view of the uncertainties in both sets of calculations, 
the agreement between the sum rule and lattice results is remarkable. For the total rate, Ball and 
Braun find 

r(B° -> p + e~v e ) = (13.5 ± 1.0 ± 1.3 ± 0.6 ± 3.6)|Kb| V" 1 . (38) 

to be compared to the lattice result in eq. ([m]) . The second error in eq. j38|) is due to the uncertainty 
in the mass of the 6-quark, the remaining errors are estimates of various uncertainties in the light- 
cone sum rule calculation and in the distribution amplitude of the p-meson. 



4 Inclusive Semileptonic Decays - V u i> 

The energy of the electron (E e ) in semileptonic i?-decays is limited kinematically to a very narrow 
window: 

m % ~ m l „ 2 3 GeV < Ee < ™b ^ 2 6 GeV (3Q) 



2m B ~ 2 

This window contains only a small fraction (about 10%) of all the b — * u decays, and it is difficult 
to make theoretical predictions for the spectrum. 

Two groups have recently jaijaposed to use the hadronic invariant mass (Mft) spectrum instead of 
the electron energy spectrum E2lE2l [| About 90% of b — > u decays satisfy Mh < Mb- The spectrum 



takes the form 



a The use of the inclusive hadronic energy spectrum for the determination of V u i, had been proposed in ref. I 
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The non-perturbative QCD effects are contained in S(y), and specifically in the parameters a% which 
can be expressed as matrix elements of composite operators: 

(H(v)\h v (iD^J(iD^ 2 )---(iD^ n )h v \H(v)) =2mBa n v^ 1 v IX2 ---v fln , (41) 

where ao — 1, a\ = and ai is given in terms of the kinetic energy of the heavy quark in the meson. 

It is proposed to determine V u b from the integral of the spectrum up to some maximum hadronic 
mass; the precision will depend on the value of the cut-off which can be attained experimentally 
(the precision on V u b is estimated to be 10-20% for values of the cut-off from Mr, down to 1.5 GeV 
or so). Much theoretical and experimental work is needed to extract the optimal results from this 
method. 



5 Power Corrections 

I now digress from the main discussion of B-decays to consider the evaluation of power corrections to 
hard scattering and decay processes. Since there are many confusing statements in the literature, and 
because the evaluation of higher order terms in the heavy quark expansion is very important in B- 
physics, it may be useful to clarify some of the key points. Although the discussion is presented in the 
context of J5-decays, it can readily be applied to other processes for which an OPE is useful and where 
the power corrections are given by the matrix elements of higher-twist or higher-dimension operators. 
The approach presented here was developed together with G. MartinelliEj, where references to the 
key papers can be found. 

Consider some physical quantity V for which the OPE allows us to write the theoretical predic- 
tion in terms of an expansion in inverse powers of mj,: 

V = C (m 6 ,//){/|O (M)K) + -4a 1 (m 6j M){/|Oi(i")N) + ■ • ■ (42) 

where the ellipsis represents higher order terms in the OPE which we will not consider further. 
For simplicity I assume here that there is only a single operator in each of the first two orders of 
the expansion; if this is not the case then there will be an additional mixing of operators which 
only requires a minor modification of the discussion below. For clarity of notation I suppress the 
dependence of the Wilson coefficient functions (d) on the strong coupling constant. 

I assume that we are interested in evaluating the 0(l/ni£) corrections in V . The theory of 
power corrections is very delicate, since they are exponentially small compared with the terms of 
the perturbation series for Co (l/m(, ~ exp(— c/ci! s (m&)), where c is a constant). Thus in order to 
evaluate the power corrections, we need to control the exponentially small tail of the perturbation 
series for Co . This raises the problems of the Borel summability and uniqueness of the sum of this 
perturbation series, and the question of renormalon ambiguities which I will not discuss further 
hcrcoE3. 

Frequently, renormalization schemes based on the dimensional regularization of ultraviolet di- 
vergences are used to define renormalized operators (e.g. the MS scheme). One consequence of 
the problems mentioned above is that higher dimensional operators (such as Oi) are not uniquely 
defined in the MS scheme; the remaining ambiguity in the matrix elements of 0\ is of 0(Aq CD ), 
i.e. of the same order as the matrix elements themselves). The exceptions to this are operators 
whose matrix elements give the leading contribution to some physical quantity, e.g. the dimension 
5 chromomagnetic operator 

A * = 3 2Mb Wh6re 4A2 = M -- M - + UJ (43) 
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and h is the field of the (static) heavy quark. 

In order to avoid the renormalon ambiguity it is necessary to introduce a hard ultraviolet cut-off 
(or subtraction scale), n, so as to provide a well-defined boundary between long- and short-distance 
contributions. In lattice calculations this occurs naturally; the lattice sparing a is such a cut-off 
(/i = a -1 ). In continuum calculations, a hard cut-off has been used in refsE30'E3. 

With a hard cut-off it is clear that the matrix elements of higher dimensional operators diverge 
as powers of fi and in themselves do not have a natural physical interpretation, e.g. 

(/KM/x)|i)~/iV (44) 

Since physical quantities cannot depend on the cut-off, the term proportional to /i™ present in the 
matrix element in eq.(^J) must be cancelled in the full prediction of V . Specifically, this term is 
cancelled by corresponding ones in the coefficient function Co, which now takes the form 

n 

Co = E + ^nE d * a *- ( 45 ) 

i=0 b i=l 

The second term on the right hand side of eq. (^5|) arises from the matching of full QCD onto 
the hqct, when the operator 0\ is included in the OPE. The important point is that the coefficient 
function Co is evaluated perturbatively and hence, in practice, this cancellation will only be achieved 
partially. In order to evaluate (or even to estimate) the 0(Aq CD /m^) corrections to V, we must 
evaluate Co to a sufficiently high order of perturbation theory to reach the corresponding precision. 
We believe that this is not the case in present calculations, and certainly has not been demonstrated 
to be so 0. 

For calculations in heavy quark physics the above discussion implies that there is no "natural" 
definition of parameters such as the binding energy A (see eq. ( f!o|) ) or the kinetic energy Ai 

X ^B) = ^(B\h(^D) 2 b\B) . (46) 

Different definitions of these parameters, due to different renormalization prescriptions for the oper- 
ators, differ by terms of O(Aqcd) and 0(Aq CD ) respectively. It is therefore of little use to compare 
values of these parameters, obtained using different definitions (often the definitions are implicit and 
must be inferred from the details of the calculations); nevertheless, this is frequently done. 



6 Non-Leptonic Inclusive Decays 

In this section I discuss two very interesting problems in the phenomenology of £?-decays, that of 
lifetimes and the semileptonic branching ratio. The discussion will use the formalism of Bigi et 
al.(see ref. E3 and references therein), developed and used by them and many other groups, in which 
inclusive quantities are expanded in inverse powers of the mass of the heavy quark, e.g. 

1927r ci [_ \ 2m^ / mf \ m b J ) 

where T is the full or partial width of a beauty hadron Hf,, 03,5 are coefficients which can be computed 
in perturbation theory and A1.2 are the parameters introduced in section |5| above. Here I will not 
rediscuss the cancellation of renormalon ambiguities present in C3, Ai and the quark mass; below 
we will consider ratios of physical quantities for which the cancellation is more transparent. An 
important feature of the general expression in eq. ( f47l) is the absence of terms of 0(l/mb), which is 
a consequence of the absence of any operators of dimension 4 which can appear in the corresponding 
OPEB. 

b In lattice calculations it may be possible to calcula*« the coefficients di to reasonably high orders by using the 
Langevin stochastic formulation of lattice field theoryE2l. 
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6.1 Beauty Lifetimes 

Using the expression in eq. for the widths one readily finds the following results for the ratios 
of lifetimes: 



1^2 = i + o(l 



! + — ) ( 48 ) 



r(A fc ) = 1 . /4(A„)-/4(2?) , /4(A b )-M 2 G (i?) 

= (0.98 ± 0.01) + (i) , (49) 



V"6 

where = — Ai and /^g = 3A2. In order to obtain the result in eq.(^), one needs to know the 
difference of the kinetic energies of the 6-quark in the baryon and meson. To leading order in the 
heavy quark expansion we have: 

2/n 3, m M B M D f M Ab -M Ac 3 M B . - M D , 1\ 
^(A b )-^(B) = — — -- — — --J . (50) 

From equation (|50|), and using the recent measurement of m AB from CDF, one finds that the right 
hand side is very small (less than about 0.01 GeV 2 ). The matrix elements of the chromomagnetic 
operator are obtained from the mass difference of the B*- and B-mesons (see eq.(fi"3|) and from the 
fact that the two valence quarks in the At, are in a spin-zero state. The theoretical predictions in 
eqs. (|8|) and ( f49| ) can be compared with the experimental measurements 

T ^* B 2 = 1.06 ±0.04 and ^|^y = 0.79 ± 0.05 . (51) 

The discrepancy between the theoretical and experimental results for the ratio t(A{,)/t(B ) in 
eqs. (|49| ) and ( pl| ) is notable. It raise the question of whether the 0(l/mf) contributions are 
surprisingly large, or whether there is a more fundamental problem. I postpone consideration of the 
latter possibility and start with a discussion of the 0(l/ml) terms. 
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Figure 8: Examples of diagrams whose imaginary parts contribute to the total rates for the decays of beauty hadrons 
(left-hand sides) and the operators they correspond to in the Operator Product Expansions. T represents a Dirac 
matrix. 

At first sight it seems strange to consider the 1 /m\ corrections to be a potential source of large 
corrections, when the 0(l/m|) terms are only about 2%. However, it is only at this order that the 
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"spectator" quark contributes, and so these contributions lead directly to differences in lifetimes for 
hadrons with different light quark constituents (consider for example the lower diagram in Fig. ^, 
for which, using the short-distance expansion, one obtains operators of dimension 6). Moreover, 
the coefficient functions of these operators are relatively large, which may be attributed to the fact 
that the lower diagram in Fig. [s] is a one-loop graph, whereas the corresponding diagrams for the 
leading contributions are two-loop graphs (see, for example, the upper diagram of Fig. ||). The 
corresponding phase-space enhancement factor is 167r 2 or so. We will therefore only consider the 



contributions from the corresponding four-quar v "™™*«~ non-wt^o- ^ti-,«,- dIa /™3 
which do not have the phase space enhancement 
these [J 



operators, neglecting other 0(l/mf) corrections 
For each light-quark flavour q, there are four of 



1 = 6 7ai (1- 7 5 )< Z q 7 "(l - 7 5 )6 ; 2 = 6(1 + j 5 )q q(l + 7 5 )6 (52) 

Ti = 6 7M (l- 7 5 )T a gg 7 ^(l- 7 5 )T a 6 ; T 2 = 6(1 + 7 5 )T a q q(l + 7 5 )T a 6 (53) 

where T a are the generators of colour SU(3). Thus we need to evaluate the matrix elements of these 
four operators. 

For mesons, following ref. o, I introduce the parametrization 

{B\O t \B)\, =mb = BJ B M B ; (BlTilB)]^ = eJ B M B , (54) 

where /i is-Jthe renormalization scale. We have chosen to use m b as the renormalization scale. 
Bigi et al. E2I prefer to use a typical hadronic scale, and estimate the matrix elements using a 
factorization hypothesis at this low scale. Operators renormalized at different scales can bc-xelated 
using renormalization group equations in the hqet (sometimes called hybrid renormalizationEj) . For 
example, if we assume that factorization holds at a low scale /j such that a s (/x 2 ) = 1/2, then, using the 
(leading order) renormalization group equations, one finds Bi = B 2 = 1.01 and e x = e 2 = —0.05 |j. 
In the limit of a large number of colours N c , Bi = 0(N®) whereas = 0(1/N C ). 
For the A b , heavy quark symmetry implies that 

(A b |0 2 |A b ) = -i^lOilAb) and (A b |T 2 |A b ) = -i(Ai,|Ti|Af,) , (55) 

so that there are only two parameters. It is convenient to replace the operator Ti, by 0\ defined by 

6 X = b%(l - 7 V ?V(1 ~ 1 5 W , (56) 

where i,j are colour labels, and to express physical quantities in terms of the two parameters B and 
r defined by 

(A 6 |6i|A,,) M=mb = -BiAblO^^m, (57) 



21/ JA b | 0l |A b)ra = _^£r (58) 

We do not know the values of these parameters. In quark models 5=1, and r = 0.2-0.5. Using 
experimental values of the hyperfinc splittings and quark models, it has been suggested that r may 
be largera, e.g. 

^!it^ = °- 9±0 - 1 ^§ 17 £^ = i. 8 ±o.5. ( 59 ) 

FT e £ 

C I use the notation of ref. L3- 

d By factorization we mean that if the Bi's and ti's had been defined at this scale (instead of mi,) they would have 
been 1 and respectively. 
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Figure 9: Theoretical Prediction of the semileptonic branching ratio and charm counting. The data points are the 
experimental results from high-energy (LEP) and low energy (i.e. at the T(4S) from CLEO) experiments. 



The lifetime ratios can now be written in terms of the six parameters -Bi,2, ei.2, B and r (as well 
as / B ): 

= J ! . I {0.02^1 + 0.0052 - 0.70ei + 0.20e 2 } 



t(B- 



t{B°) 
t(B°) 



0.98 



200 McV 
1b 



200 MeV 
+ (-0.01-0.025)?-} , 



{-O.OOBi + O.OOBa - 0.17ei + 0.20e 2 



(60) 



(61) 



where the effective weak Lagrangian has been renormalized at /1 = m&. The central question is 
whether it is possible, with "reasonable" values of the parameters, to obtain agreement with the 



experimental numbers in eq. (51). At this stage in our knowledge, the answer depends somewhat on 

up< 

Sb 



what is meant by reasonable. For example, Neubert, guided by the arguments outlined above, has 
considered these ratios by varying the parameters in the following rangesl! 



Bt,Be 



'2 


4~ 


; ei € 


1 


1" 


3' 


3_ 


3' 


3 



€ [0.25,2.5] 



200 McV 



e [0.8,1.2] 



(62) 



He concludes that, within these ranges, it is just possible to obtain agreement at the two standard 
deviation level for large values of r (r > 1.2) and negative values of 62. Lattice studies of the 
corresponding matrix elements are underway; a recent QCD sum-rule calculation has found a small 
value of r, r ~ 0.1-0.3E3. 

If the lattice calculations confirm that the parameter r is small, or find that the other parameters 
are not in the appropriate ranges, then we have a breakdown of our understanding. If no explanation 
can be found within the standard formulation, then we will be forced to take seriously the possible 
breakdown of local duality. This is beginning to be studied in toy field theories OH. 

6.2 The "Baffling" Semileptonic Branching Ratio 

This was the name given by Blok et al. El to the observation that the experimental value of the 
semileptonic branching ratio 

T(B Xev) 



B 



SL 



Ei res -► xid) + r had 
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(63) 



appeared to be lower than expected theoretically. In eq. (|62j) the sum is over the three species of 
lepton, and Th ac i and T rare are the widths of the hadronic and rare decays respectively. Bigi et al. 
concluded that a branching ratio of less than 12.5% cannot be accommodated by theoryl£j. Since 
then Bagan et al. have completed the calculation of the 0(a s ) corrections- and in particular of the 
b — > ccs component (including the effects of the mass of the charm quark)E3; these have the effect of 
decreasing B$l- With M. Neubert, we used this input to reevaluate the branching ratio and charm 
counting (n c , the average number of charmed particles per _B-decay)t3 finding, e.g. 

B SL = 12.0 ±1.0% (fi = m b ) n c = 1.20 =F 0.06 (fi = m b ) 

B SL = 10.9 ±1.0% (fi = m b /2) n c = 1.21 T 0.06 (/z = m b /2) . [ ' 

fi is the renormalization scale and the dependence on this scale is a reflection of our ignorance of 
higher order perturbative corrections. The experimental situation is somewhat confused, see Fig. ^. 
In his compilation at the ICHEP conference last yearlTJ, Richman found that the semileptonic 
branching ratio obtained from B-mesons from the T(45) is g: 

Bsl(B) = (10.23 ± 0.39)% , (65) 

whereas that from LEP is: 

B S L(b) = (10.95 ± 0.32)% . (66) 

The label b for the LEP measurement indicates that the decays from beauty hadrons other than the 
B-meson are included. Using the measured fractions of the different hadrons and their lifetimes, 
and assuming that the semileptonic widths of all the beauty hadrons are the same, one finds: 

B SL {b) = (10.95 ± 0.32)% =*> Bsl(B) = (11.23 ± 0.34)% , (67) 

amplifying the discrepancy. It is very difficult to understand such a discrepancy theoretically, 
since the theoretical calculation only involves T$l (and not T^ad for which the uncertainties are 
much larger). In view of the experimental discrepancy, I consider the problem of the lifetime ratio 
T(Ab)/r(B°), described in subsection O] above, to be the more significant one. 



7 Exclusive Nonleptonic Decays 



In this section I consider two-body nonleptonic decays of B-mesons n for which a large amount of 
data is becoming available, particularly from the CLEO collaboration!!! This is an exciting new field 
of investigation, which will undoubtedly teach us much about subtle aspects of the standard model. 
Unfortunately, at our present level of understanding we are not able to compute the amplitudes 
from first principles, and are forced to make assumptions about the non-perturbative QCD effects; 
frequently these assumptions concern factorization. These assumptions may well be wrong. Thus 
the analyses are limited to a semi-quantitative level. In this talk I will briefly describe some recent 
attempts to understand nonleptonic exclusive decays; at this stage it is not possible to endorse these 
approaches with any confidence. 



7.1 Generalized Factorization Hypothesis 

Neubert and Stech suggest an_approach based on keeping the leading order terms in the limit of a 
large number of colours (iV c )L3. Consider B — » Dn decays, for which the effective Hamiltonian is 
given by: 

Heff = {cijft) (du) (cb) + c 2 (fi) (cu) (db)} , (68) 

e Note that the rapporjeiir at the 1997 EPS conference argued that the branching ratio had been overestimated by 
the LEP collaboration 
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where, at the renormalization scale /i = mj, ci(m{,) = 1.13 and 02(7715) = —0.29, and the V — A 
structure of the current is implied, e.g. (du) (cb) = (cfy^l — r y 5 )u) (c"y M (l — 7 5 )fr)- 

For the class-1 decay B° — > D + ir~ , using colour and spinor Fierz identities we can write: 



iB°^D+7r- 



Cl 



C2 



(D+7T- I (du) (&)|B°> + 2c 2 (L>+7T- I (dT a u) {cT a b)\B } 



Following ref. we write 



^ (BD ^ = {w-\(du)\0){D+\(cb)\B°) 



and 



BD)m 



with 



a x = ci 



eg 



(BD)ir 



(69) 



(70) 



(71) 



So far this is only a parametrization of the amplitude. The two factors on the right-hand side of 
eq. ( |70| ) are given in terms of the decay constant /„■ and the form-factors for semileptonic B — > D 
decays respectively. The naive factorization hypothesis would imply that a± = c± + c 2 /N c . Neubert 
and Stech argue that the large- N c expansion may be a better guide than factorization, in which case 
we have oi = ci + 0(1/N*). 



(72) 



For the class-2 process B° — > D°7r , the relations corresponding to eqs. (69) to ( |7l|) are: 



A R0 



Ft 



(Btt)D 

a-2 



T o = (^2 + jf) (D°Tr a I (cu) (d6)|-B°) + 2 Cl (y3 7r | (cT a u) (dT Q 6)|B°) 
(D°|(cu)|0)( 7 r |(d6)|S°) 

02-?"(B7r)D 



C2 



(Btt)D 



cie 



(B-rr)D 



(73) 
(74) 

(75) 



In this case naive factorization would imply that 02 = c 2 + Ci/N c , whereas using the large iVc 
expansion a-i = C2 + Ci (1/N C + e[ B7T)D ) + 0(1/7V C 3 ), (the leading terms in a 2 are all of 0(l/N c )). 
This approximation preserves the correct renormalization group behaviour (up to corrections of 
0(1/N?)). 

Neubert and Stech propose a "generalized" factorization hypothesis for processes in which a 
large amount of energy is released, based on the large Nc expansion described above, and on the 



ancept of colour transparency applied also to class-2 decays O (for more formal arguments see also 
)0 

ai = Ci(m b ) and a 2 = c 2 (m b ) + Cci(m b ) , (76) 

where £ is a process independent parameter (for two-body decays). They study a wide variety of 
processes and conclude that a\ >2 are process independent within the available precision: 



ai = 1.10 ± 0.07 ± 0.17 and a 2 = 0.21 ± 0.01 ± 0.04 



(77) 



It is now important, not only to extend the theoretical and experimental studies to improve the 
precision and determine the range of validity of the hypothesis, but above all to try and understand 
the theoretical foundations (if any) for it. 

f Neubert and Stech alsorapply factorization assumptions to processes in which the energy of the outgoing particles 
are not necessarily large L3. 
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7.2 B -Decays to Two Light Hadrons 

The recent data from the CLEO collaboration on B-decays into two light mesons has stimulated 
many theoretical papers (see e.g. refs. E30 and for decays (inclusive as well as exclusive) with 
an rf in the final state in particular see refs.Ea). Many of the decays have suppressed tree level 
contributions, so that loop effects, which are sensitive to the presence of new physics, are important. 
As always, the principal difficulty in drawing quantitative conclusions from the experimental data 
is our inability to control the non-perturbative QCD effects. For example, control of the penguin 
contribution to the decay B — > 7r + 7r~, which may well be significant, is needed in the forthcoming 
studies of CP-violation and in the determination of the angle a of the unitarity triangle. 

Among the recent phenomenological studies of the CLEO data is a standard analysis by Ali and 
GreubL3 who include the next-to- leading order coefficient functions, and estimates of the hadronic 
matrix elements based on factorization. They find that all the B — > Kt: and tttt rates or upper limits 
can be accommodated, within this picture. Their estimates of the effects of charm quarks in penguin 
diagrams are based on perturbation theory. This has been criticized by Ciuchini et al. who stress 
that the relevanLxegions of phase space are close to the charm threshold and hence intrinsically 
non-perturbativeLl Their book-keeping of penguin effects within the OPE is briefly outlined in the 
next subsection. Concerning the phenomenological conclusions, Ciuchini et al. question whether 
the consistency of theoretical predictions and experimental measurements or bounds found in ref. l£I 
can be achieved for all processes with a single set of non-perturbative parameters. For example, 
Ciuchini et al. find it difficult to satisfy simultaneously the experimental branching ratio for the 
process B — > r)'K and the bound on the branching ratio for Bd — > 7r + 7r _ . Such debates in the 
community, on what is a very new field of study, will lead to considerable insights into subtle 
features of the Standard Model. 



7.3 Charming Penguins 

The analysis of non-leptonic 6-decays generally starts with a classification of the relevant operators 
and a calculation of their coefficient functions in perturbation theory. For example: 

Ql = (bd) V - A (uu)v- A Q\ = (bd) V -A(cc)v-A , , 

Q1 = (bu) V -A(ud)v-A Q2 = (bc) v ^ A {cd) V -A , 

plus QCD and electroweak penguin operators. The coefficients of Q\'^ are of O(l) whereas those of 
the remaining operators are of 0{a s ) and are generally small. i_. 

Ciuchini et al.113 have recently emphasized the point, previously made by Buras and Fleischer l2I 
that the matrix elements of the operators Q\'^ have penguin-like contractions (when the up or charm 
quark fields are Wick contracted). For the case of the charm quark, the effects of these "charming" 
penguins are enhanced in decays where the emission diagrams are Cabibbo suppressed relative to 
the penguin diagram. Although the precise definition of charming penguins (i.e. the separation of 
penguin effects between the explicit contributions of the further operators Q3,... and those in the 
matrix elements of the operators Qi^) is a matter of convention, they do contain non-perturbative 
QCD effects which must be included in the calculations of decay amplitudes. From their recent 
analysis, Ciuchini et al. conclude that large contributions are likely from charming penguins, e.g. 
that it is necessary to include them in order to obtain branching ratios of order 1 s^_L0 5 for the decays 
B + — ► K°n + and Bd — > K + n~ as recently found by the CLEO collaboration EJ. By using recent 
results from CLEO, Ciuchini et al. constrain many of the relevant hadronic matrix elementa_and are 
able to make predictions for a large number of processes which have not yet been observedLO. Some 
of these processes, including B — > pK, are predicted to have branching ratios close to the current 
experimental bounds, and will provide an excellent testing ground for our understanding of penguin 
effects. 
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7.4 Lattice Calculations 

The amplitudes for non-lcptonic decays of hadrons {H — ► h\h 2 ) are difficult to evaluate in lattice 
simulations in principle. Maiani ancLTesta pointed out that in lattice simulations one obtains the 
following average of matrix elements E3: 

\{ in {hih 2 \H eff \H) + out (h!h 2 \H eff \H) ) , (79) 

and hence no information about the phase of the final state interactions. They also show that from 
the large time behaviour of the correlation functions one can extract the unphysical matrix element: 

(MP/u - 0)h 2 {Ph 2 = 0)\H ef f\H) . (80) 

Together with chiral perturbation theory, this may be useful to obtain good estimates of the ampli- 
tudes for kaon decays, but it isjaot very useful for £?-decays. 

The Maiani- Testa theoremLa implies that it is not possible to obtain the phase of the final state 
interactions without some assumptions about the amplitudes. The importance of developing reliable 
quantitative techniques for the evaluation of non-perturbative QCD effects in non-leptonic decays 
cannot be overstated, and so attempts to introduce "reasonable" assumptions to enable calculations 
to be performed (and compared with experimental data) are needed urgently. Ciuchini et al.EaS 
have recently shown that by making a "smoothness" hypothesis about the decay amplitudes it is 
possible to extract information about the phase of two-body non-leptonic amplitudes. Studies to 
sec whether their proposals are practicable and consistent are currently beginning. 

8 Conclusions 

During recent years, through the combined work of experimentalists and theorists, there has been 
enormous progress in our understanding of heavy quark physics. In this talk I have reviewed some 
of this work, and underlined a few of the areas in which further progress is urgently needed in 
developing control of the non-perturbative QCD effects. Among the main points were: 

• The lattice community will continue to refine the computations of a wide variety of physical 
quantities (e.g. leptonic decay constants, B-parameter of B-B mixing, form factors of semilep- 
tonic decays, the parameters of the hqet), reducing the systematic errors, in particular those 
due to the use of the quenched approximation. 

• V u b can be determined by studying-the hadronic invariant mass spectrum in inclusive semilcp- 
tonic decays. The authors, of ref. EJ estimate that this will lead to an error of about 10% in 
V u b whereas those of ref. EiJ that the error will not exceed 10-20%. 

• We need to understand whether the discrepancy of the theoretical prediction for the ratio 
r(Ab)/r(B) with the corresponding experimental measurement is due to some effect which 
can be controlled (such as a matrix elements which is larger than currently expected) or to 
violations of local duality. 

• We need progress in understanding whether power corrections to hard scattering and decay 
processes can be controlled numerically. 

• New theoretical ideas are urgently required to interpret quantitatively the wealth of data on 
non-leptonic exclusive decays from current and future experiments. In particular progress 
in this area is needed for studies of CP-violation at the forthcoming 6-factories, and in the 
attempts to fix the angles of the uniitarity triangle. 
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